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Abstract. We consider an interacting particle system on the one dimensional 
lattice Z modeling combustion. The process depends on two integer parameters 
2 < a < M < 00. Particles move independently as continuous time simple 
symmetric random walks except that 1. When a particle jumps to a site which 
has not been previously visited by any particle, it branches into a particles; 2. 
When a particle jumps to a site with M particles, it is annihilated. We start 
from a configuration where all sites to the left of the origin have been previously 
visited and study the law of large numbers and central limit theorem for rt, the 
rightmost visited site at time t. The proofs are based on the construction of a 
renewal structure leading to a definition of regeneration times for which good tail 
estimates can be performed. 

RESUME. On considere un systeme de particules en interaction sur Z mod- 
elisant les particules incandescentes d'un mecanisme de combustion. Le processus 
depend de deux parametres entiers 2 < a < Af < cxd. Les particules se deplacent 
independamment selon des promenades aleatoires simples symetriques a temps 
continu, mises a part les interactions suivantes: 1- quand une particule saute vers 
un site qui n'a jamais encore ete visite, elle branche et fait place a a particules; 2- 
quand une particule saute vers un site abritant M particules, elle disparait. On 
demarre d'une configuration ou seuls les sites a gauche de I'origine ont deja ete 
visites et on etudie la loi des grands nombres et le theoreme de la limite centrale 
pour rt, la position du site le plus a droite visite a I'instant t. Les preuves re- 
posent sur la construction d'une structure de renouvellement associee a des temps 
de regeneration dont les queues peuvent etre convenablement estimees. 



1. Introduction 

The method of regeneration times has been very successfully applied to problems 
of random walk in random environment (for example see |17| . |18| . j^). In a one- 
dimensional setting it was already used via a renewal structure in jH]. On the other 
hand, a concept similar to that of regeneration times, known as duster structure, 
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has been well know for many years in the context of mechanically interacting one- 
dimensional dynamical systems of particles (also see pH] and references therein). 
In this article, we extend the method of regeneration times to an interacting particle 
system. The system we consider is a one dimensional stochastic model of combustion. 
Heat particles move by symmetric nearest neighbor random walks on sites {x G 
Z : rr < r} of the integer lattice, with r representing the position of a flame front. 
To the right of the flame front is a propellant. The first heat particle which reaches 
the propellant at r + 1 immediately branches into a > 2 particles, and the front 
moves one step to the right. We also include an upper bound M on the number of 
particles at each site, so that if a particle tries to jump to a site with M particles, 
it is immediately killed. We show that the front in the model moves ballistically to 
the right, or, more precisely, prove a law of large numbers for r^. The next question 
one might ask is of the fluctuations about the law of large numbers. We prove here 
that these are Gaussian, with a central limit theorem for t~^l'^{rt — vt). 

For M = oo a discrete time version of the system we are considering has appeared 
in the literature under the enigmatic name "frog model", and laws of large numbers 
for the position of the front were proved, also in higher dimensions, using methods 
based on the sub-additive ergodic theorem [2], [H], @], ^21, ^H]- The continuous 
time case under the name stochastic combustion process was treated in ^3], [l5j . 
Note that the M = oo case is sub-additive, but the M < oo case is not. Our 
main purpose here is to develop new methods to study such models, and we are 
especially interested in the fluctuations of the fronts. If nt represents the number 
of particles at the front rt at time t, then rt moves to the right at rate nt, and 
hence rt = /q risds + Mt where Mt is a martingale. The standard approach to the 
law of large numbers would then be to show that, as observed from the front, the 
system has an ergodic invariant measure /i, and j^risds — > E^[n]. The central 
limit theorem would be proved by showing that the time correlations of rj{s) decay 
fast enough. However, very few methods exist for proving uniqueness, or ergodicity, 
of invariant measures of such systems. So such approaches appear to have limited 
applicability. Furthermore, although it seems intuitive that the correlations of r/(s) 
decay quickly-likely exponentially fast- it is not at all apparent how to prove it. In a 
future article, we will consider the case M = oo which can be handled by essentially 
the methods developed here, though with major technical modifications. 

The combustion process that we consider is related to deterministic reaction- 
diffusion equations of the form ^ = + f{u). The two key differences are the 
discreteness of the variable-we have a number of particles as opposed to a continuum 
variable tt-and the stochasticity. The effect of discreteness and/or stochasticity on 
the traveling waves of reaction-diffusion equations is a question that has not received 
the attention it deserves, as these effects are likely present in real systems. Therefore, 
we believe it is crucial to develop methods to study such systems. In the literature 
of reaction-diffusion equations one separates several cases according to the behavior 
of / near zero. An / which vanishes on [0, 9] for some > with f[u) > for u > 9 
is said to have a combustion nonlinearity with ignition temperature cutoff 9 (see, 
for example Note that the discreteness of the particle models of the type we 
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consider makes them essentially of the combustion type with an ignition temperature 
of one particle. 

In the second section of this article, we define the combustion process and state the 
main result (Theorem^. In Section 3, we define an auxiliary and a labelled process, 
which will be needed to define the renewal structure leading to the regeneration times. 
The labelled process can be understood as a combustion process where particles are 
labelled so that if at a given site a particle has to be killed, it is the one with the 
smallest label. The auxiliary process moves ballistically to the right and is coupled 
to the labelled combustion process in such a way that it is always to the left of 
the right-most visited site r^. In Sectional these processes are used to define the 
renewal structure, and the corresponding regeneration times. Then in Section it is 
proved that the regeneration times have finite moments of order 2, under appropriate 
assumption on the threshold M. In Section El we complete the proof of Theorem ^ 



We define a stochastic process describing the dynamics of particles on the lattice 
Z which move as rate 2 continuous time simple symmetric random walks and branch 
and are killed at rates depending on the configuration of neighboring particles. The 
branching and killing depend on natural number parameters 2 < a < M. Each 
particle performs a continuous time symmetric simple random walk independent of 
the others, with two twists: There is a position r S Z which is the rightmost visited 
site. When a particle jumps right from this site, it branches into a particles, with the 
result that there are a particles at the new rightmost visited site, r + 1. In addition, 
we allow at most M particles at any site, and maintain this requirement by killing 
any particle that attempts to jump to a site with M particles. 

The state of our system is 



2. Combustion process 



n := {(r,r/) 
The infinitesimal generator is 

Lf{r,r]) = 



rGZ,r?G{0,...,M}^-''^-^''^>}. 



^ r]{x){f{r,T^yr]) - f{r,r])) 



x<r,y<.r,\x—y\=l 




(1) 



(2) 



8 < M < oo. 



(3) 



We can now state the main results. Throughout the sequel we will use the notation 
7]{t, x) to denote the number of particles at time t > at site x in the stochastic 
combustion process. 
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Theorem 1. Suppose the process is started with initial conditions r = 0, r]{0,0) £ 
{1, . . . , M} and 77(0, x) G {0, . . . , Af } arbitrary for x < 0. 

(i) (Law of large numbers). There exists v,Q <v < 00, which does not depend on 
rj, such that a.s., 

lim rt/t = V. (4) 

f— »oo 

(a) (Central limit theorem). There exists a, < a < 00, which does not depend on 
the initial condition {r]{0,x) : x G Z}, such that 

e'/' {r,-H - e'^vt) , t>0, (5) 

converges in law as e — > to Brownian motion with variance < o"^ < 00. 

To prove this theorem we will define a renewal structure for the right-most visited 
site Tf, in terms of regeneration times, and then will show that such times have finite 
second moments. 



3. Auxiliary and labelled processes 

In this section, we will define an auxiliary process and a labelled process, both 
coupled to the combustion process. The auxiliary process {r^ : t > 0}, will take values 
on Z and will have two fundamental properties: under certain initial conditions on 
the combustion process such that ro = ro, it will always be to the left of the right- 
most visited site {rt : t > 0} and its dynamics is independent of the particles to 
the left of ro in the combustion process if the appropriate rule for killing is defined. 
These properties will help us define the renewal structure for the right-most visited 
site {rt t > 0} of the combustion process and will give us easily a lower bound for 
the limiting speed v of rt- This is the content of Lemma ^ The labelled process, 
corresponds to a process with a state space larger than the combustion process, where 
an explicit rule for killing is given: particles are labelled, and every time a particle 
has to be killed, it is the one with the smallest label which is killed. 

Let us first introduce some general notation. Define ^ := {1, 2, . . . , a}, B ■.= 'LxA 
and consider a set of independent continuous time simple symmetric random walks 
{Yx^i : G B}, each one of total jump rate 1, and such that Yx^i{Q) = x for each 
i ?i A. We want to endow the set of indices B with the lexicographic order. 

Consider a set of independent continuous time rate 2 simple symmetric random 
walks Yx^i, with label (x, i), where x G Z and i G {1, . . . , a — 1}. Each random walk 
Yx^i starts at site x for every z G {1, . . . , o — 1}. We order the labels by 

(x^i) < {x' ,i') if X < x' or x = x' and i < i' . (6) 

3.1. Auxiliary process. Let r G Z and for each z G Z, define B^ as the set of M 
labels not exceeding (z,a — 1) in the order (EJ. Define Ar^z to be the labels in Bz 
with X > r. 

We define a sequence of waiting times. Let fo := and define ui as the first 
time one of the random walks {Yz^i : {z,i) G ^r,r}, hits the site r + 1. Next, define 
V2 as the first time one of the random walks {Yz^i : iz,i) G Ar^r+i}, hits the site 
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r + 2. In general, for A; > 3, we define fk as the first time one of the random walks 
{Yz^i : {z,i) £ Ar^r+k-i}, visits the site r + k. Finally, for n E N, let 

f[:=r + n, if Ek=o''k < t < Ektl^^k, (7) 
where the superscript r in f[ indicates that fg = r. 

Lemma 1. There exists an a > which does not depend on r such that a.s. 

lim fl/t = a. (8) 

t^oo 

Proof. Once one notes that for every j > the random variables {yu+j '■ k>\} are 
independent whenever I > M/a, the proof is a simple exercise. □ 



3.2. Labeled process. We enlarge the state space of the stochastic combution pro- 
cess so that particles carry labels which tell us where they originated. 

We will want to keep track of where particles came from, even after restarting at 
stopping times. Hence each particle will have a starting position z € Z and label 
(x, i), X G Z, i G {1, . . . , a— 1} describing its birthplace. We will allow the possibility 
that z ^ X. 

At time 0, we have an r G Z representing the rightmost visited site, and a subset 
X(0) of the labels (x, i) with x < r, representing the set of labels of live particles 
at time 0. To each one of these labels is assigned a position z = Zx^i{0) < r 
which is the position at time i = of that particle. The position at time t is 
Zx,i{t) = Yx,i{t) + z- X. 

To keep track of the killing in this process, we make a rule that whenever a particle 
jumps to a site with M particles, the particle at that site with the smallest label is 
killed and the corresponding label is removed from the set T of labels of live particles. 

When a particle jumps to site r + 1, the labels {(r + 1, 1), . . . , (r + 1, a — 1)} with 
corresponding to particles with initial positions r + 1, are added to the set of labels 
of live particles. The time this happens will be denoted pi. These particles then 
have trajectories which is equal to Yr+i^i{t — pi) foi t > pi. 

Similarly, for > 2, + • • • + will be the first time a particle jumps to r + A; 
and at that time {(r + /c, 1), . . . , (r + /c, a — 1)} are added to T, with trajectories 
Zr+k,i{t) = Yr+k,iit — Pi — ■ ■ ■ — Pk) for t > pi + ■ ■ ■ + Pk Until such time as the 
particles are killed and their label removed from the set of labels of live particles. 

We denote by X{t) the set of labels of live particles at time t and by Z{t) = 
{Z[x,i){t) ■ i^ji) £ the positions of the corresponding random walks. 

To avoid pathologies it is useful to insist that initially the set of labels of live 
particles includes at least one with x = r. The rightmost visited site r is the supre- 
mum of the x over the collection of labels, rt = sup{x : (x,i) G 1{t)}. Let us 
formalize the above discussion. Call L the triples {r,T,Z) formed by an integer 
r G Z, a set of labels I C {{x,i) : x < r,l < i < a — 1} and position function 
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Z : I —> {. . . ,r — 2,r — l,r} taking values in the integers smaller than or equal to r. 
We now define 



(r,Z,Z) C L : max x = r, max > '^z(xi)i^) ^ ^'^ 



which will be the state space of our process. 

The labeled process starting from w = (r,X(0), Z(0)), is now defined as the triple 
Wt = {{rt,I{t), Z{t)) : t > 0}, defining a strong Markov process taking values on S, 
and with a law given by a probability measure defined on the Skorohod space 
^([0,oo);S). 

The combustion process described in section 1 is the particle count r]{t) := {r]{y, t) : 
y < rt} where r]{y,t) = Yl(xi)GJ{t) ^{^x,i{t) = u)- We will occasionally use the more 
explicit notation notation rj^i^t) and rj^{y,t) instead of ri{t) and rj{y,t) respectively, 
with the understanding that 77«,(0) is the particle count of the initial condition w. 

We already defined pi to be the first time that one of these particles hits r + 1. 

Lemma 2. Suppose that (r, 1), . . . , (r, a — 1) G 2^(0), all initially at r. Then pi < vi. 

Proof. The dynamics of the a particles which start at r is the same in the labeled 
process as in the ones we look at in the auxiliary process. Indeed, these are the only 
particles considered in the auxiliary process up to time ui, while the labeled process 
may have many others, each of which has a chance to be the first to hit r + 1. □ 



Recall that Pk + • • • -\- Pi is the first time one of the labeled particles in the labeled 
process hits r + k. Note of course that the hitting could be done by one of the a new 
particles created at r + i at time pi, i < k. 

Lemma 3. Suppose that (r, 1), . . . , (r, a — 1) G 1(0), all initially at r. Then pk < I'k- 

Proof. Before giving a general proof, we describe the special case ofM = a = k = 2 
where the idea is more transparent. 

Case M = a = k = 2: Note first of all that if a = 2, we have i = 1 always 
and hence we can drop the i in the labels. The labeled process starts with one 
particle with label r at r, possibly one other particle at site r, with a label smaller 
than r, and other particles with labels smaller than r at arbitrary positions to the 
left of r. At time pi we have one particle labeled r + 1 at r + 1, which, after that 
time, has trajectory Zj.+i{t) = Yr+i{t — Pi), and another particle, labeled r at some 
site X < r + 1, which, after that time, has trajectory Zr{t) = Yr{t). Note that 
neither particle can be killed before time pi + p2 because they are the two with 
the highest labels until that time. There could also be other particles with other 
labels at positions x < r. Denote by r^+i = inf{t > pi : Zr+i{t) = r + 2}, by 
Tr = inf{t > Pi : Zr{t) = r + 2} and by Tothers the first time one of the others hits. 

Then P2 = m:m.{Tr+l,Tr, Tothers] - Pi < m.m.{Tr+l,Tr} - Pl. 
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On the other hand, 1^2 = mm{(Tr, (Jr+i} where 0"^ = inf{t > : Yj.{t) 
(Jr+i = inf{t > : Yr+i{t) = r + 2}. Now 

Tr+i = inf{i > pi : Zr+i{t) = r + 2} = mi{t > pi : Yr+i{t - pi) = 
= M{t > : Yr+iit) = r + 2} + pi = cr^+i + pi. 

and 

Tr = inf{t > Pi : Zr{t) = r + 2} = inf{t > pi : Yr{t) = r + 2} 
< mf{t > : Yr{t) =r + 2}+pi=ar + pi. 

The result follows. 

General case: Let us examine the configuration in the labeled process at time 
Pi + - • • + pk-i- It has a — 1 particles at site r + k — 1 with labels {r + k — 1,1), ... ,{r + 
k — 1, a — 1) plus one additional particle with an unknown label, the one which hit 
r + k — 1. An additional M — a + 1 particles (which might include the additional 
particle of unknown label which hit r + A; — 1) with the next highest labels after the 
first a — 1 have some unknown positions, and there may in addition be any number 
of other particles in the configuration.. In the time interval [0, pi + • • • + Pk~i] none 
of the M — a + 1 particles has hit r + k. Their trajectories in the time interval 

[pi-\ \-Pk-i,Pi-\ \-Pk] areY^j.+k-2,a-i)it- Pi-\ \-pk-2),---- The definition 

of ffc, on the other hand, involves the particles with these M leading indices, but 
without the time shift. Hence the first time that any of the first a — 1 particles 
hits is identical to that of the auxiliary process, the first time that any of the next 
M — o + 1 particles hits is greater in the auxiliary process, and one of the other 
possible particles in the labeled process could be the one which hit first, making the 
time shorter still. □ 

3.3. Enlarged process. We have seen that the labeled process is defined in terms 
of a set of labels I(t), for t > 0, with the property rt = sup{x : {x,i) £ lit)} < 00. 
Whenever a particle with a label in I{t) is killed, this label is removed. In the 
enlarged process we keep track of the killed particles as well. Define for each t > 0, 
the set of labels of all activated particles up to time t in the labeled process, 

i{t) = uo<s<a(.s). 

Consider the set Z{t) = {Zx^i{t) : {x,i) £ ^(t)} of all the corresponding random 
walks at time t. The enlarged process is defined as the triple wt = {(rj,X(t), ^(t)) : 
t > 0}. 

Consider now the set of labels TZ{t), obtained after removing from I{t) all labels 
(x, i) with X < r = sup{y : {{y, i) £ 2^(0)}. We define for y < rt the particle count 

ay,t)= l(%,o(i) = y)- (10) 

Also let C{t) be the set of labels obtained after removing from T{t) all labels (x, i) 
with positions x > r. We define for y < rt the particle count 

Hy,t)= KZi.,^){t) = y). (11) 

{x,i)£C{t) 



= r + 2} and 
r + 2} 

(9) 
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We similarly define C{t) as the set of labels obtained after removing from T{t) all 
labels {x,i) with positions x >r, and the corresponding particle count for y < n as 

^{y,t)= Yl = y). (12) 

(x,i)G£(t) 

Lemma 4. (i) For every initial condition w €zS and every t > and x < rt 

(t>{x,t) <^{x,t). (13) 

(ii) For every w = (r,T(0), ^(0)) € S with labels (r, 1), . . . , (r, a — 1) at site r and 
corresponding initial positions j)(0) = r for I < i < a — 1, 

r[ < rt. (14) 

(Hi) For every w = {r,I(0), Z{0)) G S, the processes (f) and fj" are independent. 

Proof. Part (i) follows directly from the definitions. Part (ii) is a consequence of 
Lemma 131 Part (iii) follows from the observation that (f) and r[ are defined in terms 
of the random walks {Y^^i : x < r} and {Y^^i : x > r} respectively, which are 
independent. □ 

4. The renewal structure 

Consider the enlarged process wt with its natural filtration Tt with an initial 
condition wq having particles with labels (r, 1), . . . , (r, a — 1) at site r, and any 
allowable configuration of particles with labels to the left of r. Let a = limt_>oo h/t 
be as in JH)! and choose any < a' < a. Define 

U := inf{t > : fj - r < [a't\ } 

where [x\ denotes the greatest integer less than or equal to x. 

Define V as the first time one of the particles initially to the left of r — 1 hits 
{x : X > [a't\ + r}, 

y := inf{t > : sup ^^^(a;, t) > 0}. 

x>[a't]+r 

U and V are stopping times with respect to {Tt ■ t >0}. Let 

D ■.= mm{U,V}. (15) 

For each y G Z, let 

Ty := inf{t > : rt > y} 

denote the hitting time of y by the rightmost visited site in the labeled process wt. 
We will also need the first time U and V happen after time s > 0, 

Uoes:= inf{t > : - r, < [a't\}, (16) 
y o := inf{t > : sup 4>^^{x,t) > 0}, (17) 

x>[a't]+rs 

and D o Og := min{^7 o Og^V o 0^]. 



STOCHASTIC COMBUSTION MODEL 



9 



Choose an integer L := M and define sequences {S^ : A; > 0} and {D^ : A; > 1} of 
^j-stopping times as follows. Start with 5o = and Rq = r. Then define 

5i := Tro+l Di := D o esi+ Si, Ri := roi, 

and, for > 1, 

Sk+i '■= Tr^+l Dk+i -.= 9s^j^^ + Sk+i, Rk+i '■= ^Dk+i- 

Note that these times are not necessarily finite and we make the convention that 
Too = oo. Similarly, we define Uk ■= U o 6$^^ + Sk and Vk := V o Os^ + Sk for k > 1. 
Let 

K := inf{/c > 1 : 5^; < oo, Dk = oo}, 
and define the regeneration time, 

K := Sk. (18) 

Note that k is not a stopping time with respect to 

Denote by Q the information up to time k, defined as the completion with respect 
to of the smallest cr-algebra containing all sets of the form {k < t} f\ A, A ^ 

In section we will show that IE^„[k^] < oo and hence, in particular 

< oo] = 1. (19) 



Proposition 1. Let A he a Borel subset o/ D([0, oo); Jl). Then, 

[r„,^C('^ + -)^A\g]= WasM-) <^A\U = 'x], 
where oSq denotes a configuration with a particles at and none anywhere else. 
Proof. We have to show that for any B £ Q, 

F^[B,{T^r^C{K + -)£A}]=F4B] FasM-)^^ \ U = oo]. (20) 
Now, using Hlfj|l . 

F^[B, {T^rXiK + •) G A}] = F^[{k < cx)}, B, {r_,,C(^ + •) e A}] 

oo 

= J^P^ [{Sk < oo,Dk = oo},B,r_,^C('^ + •) G A] 



k=l 

oo 



^^P^[rs, =x,Sk< oo,A = oo,S,r„^C(5fc + -) e A]. (21) 



k=l x€l 



Prom the definition of Q there is an event B^ G J^Sk such that B^ = B on k = Sk- 
Therefore, we can continue developing ipT|) to obtain, 

oo 

= X^X!^"'" [^•5fc =^^Sk < oo,Dk = oo,Bk,T-xC{Sk + ■) e A] 

k=oxez 

= '^Ey, \^lrs^=x,Sk<oo,Bk^w [Dk = 00,T_^C{Sk + ■) £ A \ f's^ 
k,x 
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where E^, is the expectation defined by P^. But on the events Sk < oo and rs^. = x, 
we have by parts (i) and {ii) of Lemma 01 that 

Cv.{Sk + -)=Va5S-) (22) 

when Uk = Vk = oo, and that "HaSsi^ (") is independent of the configuration of particles 
to the left of x. Indeed, part (i) of Lemma|3 and the event Vk = oo, imply that the 
particles with initial positions z to the left of x, never hit the line [a't\ + x. And 
part (ii) of Lemma HI and the event Uk = oo, imply that the front is always to 
the right of the same line. Hence, there is no effect of the particles initially to the 
left of X on the front r^, so that CwiSk + ") = VaS^i')- Then, (|22|l combined with 
the independence of Uk and Vk given Ts^. by part {in) of Lemma ^ the translation 
invariance, and the strong Markov property imply that on the events Sk < oo and 
rSk = X, 

Fu, [Uk = oo,Vk= oo, T-^CiSk + •) e ^ I ^s,] 
= Fu, [Uk = oo, T_^r]a5^ (•) e ^ I ^s,] F«, [Vfc = oo I fs^] 
= FaS,[U = oo,r?(-) G A]F4Vk = oo | Ts,]. 
Summarizing, we have, 

F^K < oo] F4B,T^r^C{K + -) G A] 

= FaSo[U = ^,V{-) ^ A]Y,^4Vk = ^,rs, = x,Sk < ^,Bk]. (23) 

Letting A = gives 

[k < oo]P^ [B] = Fa5o [U = oo]Y, [V^fc = oo, rs, =x,Sk<oo, Bk] . (24) 

k,x 

(1^ and ^ together imply (HH). □ 

Now define ki < K2 < • • • by ki := k and for n > 1 

K„+i := K„ + k{wk„+-)- (25) 

where k{wk„+-) is the regeneration time starting from Wk„+- and we set Kn+i = oo 
on = oo for n > 1. We will call ki the first regeneration time and Kn the n-th 
regeneration time. 

For each n > 1 the cr-algebra, Qn will be the completion with respect to F^ of the 
smallest cr-algebra containing all sets of the form {ki < ti} n • • • H {k„ < i„} n A, 
AGTt„. Clearly e?! = g. 

Lemma 5. {U = oo} G Qi. 

Proof. Note that {ki = oo} is a null event for P.^, and hence, since Qi is complete, it 
is enough to show that {U < oo} n {ki < oo} £ Qi. 

For notational convenience, we write f instead of f. * . Note that whenever 
?7 < oo, 5fc < oo and rjj > rs^ happen, for some /c > 1, then necessarily ki > Sk- 
Indeed, from the observation that ^5^+. = f^, we see that if the events U < 00, 
Sk < 00 and rjj ^ fSk happen, then we must have that Uk < 00, and hence Dk < 
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OO. By summing over the intersection with {ni = Sk} we see that it follows that 
{[/ < 00} n {ki < 00} n {fu > ^Ki} = 0. So 

{[/ < 00} n {ki < 00} = {fu < "Tki} n {ki < 00}. (26) 

Since ffj < rt implies that C/ < 00, it follows that {fu > r^-^} n {fu < rt} H {ki < t} 
is empty. Hence {fu < n} f) {ki < t} = {fu < J n {fu < n} n {ki < t} = {fu < 
Tki} n {ki < t}. Thus 

00 

{fu <r^,}n {ki < 00} = \J{fu< rn} n {ki < n}. (27) 

71=1 

The result then follows from the fact that {fu < rt} G for each t > which is a 
direct consequence of the construction of the processes. □ 

Proposition 2. Let A be a Borel subset 0/ D([0, 00); $7). Then, 
Proof. Let ip : D[[{), 00); S) -^[[0, 00); S) be the map given by, 

^{w){-)=T^r^w{K^{w{-)) + -). (28) 

Then note that, Qk+i is generated by Qi and ■0"^ (^fc), and that the cx-algebras Qi 
and ip'^iQk) are independent. The proof of this theorem now follows from the above 
observations, induction on n G Z"*" using Proposition ^ and Lemma El □ 

Proposition 3. Let u; G S. (i) Under P^, ki, K2 — ki, ^3 — K2, . . . are independent, 
and K2 — Ki, — K,2, ■ ■ ■ are identically distributed with law identical to that of ni 
under FaSollU = 00]. (ii) Under F^j, r.A«:i, r(«.j+.)^«.2 - r^^^ , r(«2+.)A^3 -r^^,... are 
independent, and r(^f^^_^_.-j;^i^^ — r^j , t^^^+Oaks ~ '^k2) ■ ■ ■ '^'^^ identically distributed with 
law identical to that of r^^ under Pa5o["l^ = ^] ■ 

Proof. This follows directly from Proposition |21 □ 

5. Expectations and variances of the regeneration times 

5.1. Estimates for the auxiliary process. In this subsection we obtain some 
estimates for the auxiliary process. The process will always start from a^o and we 
denote the corresponding measure on the trajectories of the auxiliary process by P, 
and expectations by E. 

Lemma 6. For 1 < p < M/2 and j > M, there exists a constant C = C{p) < 00 
such that 

E[\u,\P]<C. (29) 

Proof. For j > M the auxiliary process has its full complement of M particles 
attempting to hit j + 1. Uj is then the minimum of 71, . . . ,7m which are the hitting 
times of 1 of M random walks starting in {—M/a < x < 0}. Standard estimates 
for random walks give Phi > t] < Ct-^/^ for some C < 00. Hence P[vj > t] < 
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Lemma 7. For 1 < p < M/2, there exists a constant C = C{p) < oo such that for 
all initial conditions w, and all t > 



P[t <U <oo]< Ct^^P/'^\ 



(30) 



Proof. By translation invariance we can assume r = 0. Let ti be such that [a'iij = 
M. Then, when t>ti, we have that 

P[t<U <oo]<P [fl > M, Us>t {r° < [a's\ }] . (31) 

But, if < [a's\ for s > t, then X]}=i^ — ■^j which in turn imphes that, 
n Sj=i — ^ fo'^ some n > [a't\ . Similarly, ff^ > M, implies that '^jLi ^ 
(M + l)/a'. Therefore, whenever [a't\ > M + 1, we have 



P[t < [/ < oo] < P [E^il-. < ^,[Jn=la'ti [iZWj > i.} 



j=Af+l — 



(32) 



Now remark that for j > M , the random variables Vj are identically distributed and 
have finite moments of order p < M/2 by Lemma El Each has expected value 1/a. 
Define jj = vj — l/a. Choose t2 as any real number such that P = — ^ — [^'t2\ ^ ^• 
Then, whenever t > t2, if N = [a't\, we have that, 



P[t<U<^]< P sup„>^ i Ei=A/+i 7j > f3 



(33) 



Recall that for each < i < i, £ = [M/a\ + 1, the random variables {fM+i '■ k >!} 
are independent. Observe from (|33|l that, if /5' = (3a/M , 



p[t<u<^\< EtlP 



SUp„>Ar 



1 Y^"- 



7, > 



(34) 



Suppose Xi,X2, . . . are independent and identically distributed random variables 
with mean 0. By Kolmogorov's inequality, 

^[SUP„>7V ^ E"=l Xi>e]< X^^o P [sup2fcjv<n<2fe+lJV W ^^^=1 ^» > 

< Er=o(2'iVe)-^^II E ^ill (35) 

Now, for p > 2, if E[\Xi\P] < oo then ^[| Eil^^^il^] < C(2^'+^A^)p/2 for some 
C < oo (see item 16, page 60 of p^), and hence for another C < oo, 

(36) 

□ 



Applying to (EU, by Lemma El we obtain (jHUJ. 



5.2. Estimates for the enlarged process. We start with a few standard estimates 
for hitting times of random walks. 

Lemma 8. Let {Xt : t > 0} be a continuous time simple symmetric random walk on 
TL, of total jump rate 2, starting from x < —1, c > 0, and 



T := mi{t >0:Xt> [ct\}. 



(37) 
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P[t = oo] > 



Then 

'l - exp{(l + x)6lc} x<-2 

exp{-2/c}(l - exp{-6'c}) x = -I 

where 9c > is the nonzero solution of c9 — 2(cosh^ ~ 1) =0, and there exist 
< C, C" < oo such that 



P[t <T <oo]< C"exp{-Ct}( exp{-6'c|2;|/2} +exp{-t/(|j;|/(2t))}), (38 



where I(n) = 2 + nsinh — \/4 + is the rate function for a random walk. 

Proof. For 6* G M, exp{9Xt — 2{cosh0—l)t} is a martingale. By the optional stopping 
time theorem, 

E[exp{eXrAn - 2(cosh 6* - l)r A n}] = exp{6'3;}. (39) 

Now, XrAn < [cT A nj , hence OXrAn - 2(cosh 6 -1)t ^n < {c6 - 2(cosh 6 - l))r A n. 
It follows that \i 9 > 9c, we can apply the bounded convergence theorem in equality 
(|39|l taking the limit when n — > oo, to conclude that, 

E[1{t < oo) exp{9Xr - 2(cosh 9 - l)r}] = exp{6ix}. (40) 

When r < oo, we have Xj- = [crj and therefore 9Xt- — 2(cosh^ — l)r > {c9 — 
2(cosh0 — 1))t — 9. Letting 9 [ 9c in and applying the bounded convergence 
theorem, we obtain, 

P[t < oo] < exp{(l + x)6lj. 
To get the bound for x = — 1, note that the probability that the random walk does 
not move before time t = 1/c is exp{— 2/c}. Then use the Markov property and the 
result for x = —2. 
To prove 

P[t<T<oo]< P[Xt > B]+P[t<T <oo,Xt< B] (41) 

Now 

P[Xt >B]< exp{-tI{{B + |x|)/t)} 
and by the strong Markov property and translation invariance, 

P[t<T< oo,Xt <B]< P_(Letj_B)[T < oo] < exp{-([ctj - B- l)9c}. 

Choosing B = {[ct\ +x)/2, on the above inequalities, and using the convexity of the 
rate function I, gives l(38|) . □ 

Lemma 9. For each a' < a, there exists C < oo such that for t > I, and all w £ E, 

F^[t <V <oo]<Cexp{-Ct}. (42) 

Proof. Note that in the worst case in which there are M random walks at each site 
to the left of the origin, we get the bound. 



[t<V<oo]<MY^ P-,[t < r < 



ooj 

x=-l 
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where r is defined in display l(37|l . On the other hand, it is true that 
e"^^'^/^""'"^-'^ < {n + 1) X^fc^o ^"^^'^^^ This estimate, the previous inequahty 
and inequahty ((38|) of Lemma |H1 give us the result. □ 

Prom Lemmas and El we have 

Corollary 1. For each < p < M/2 there is a C < oo depending only on p, M and 
a' such that for all initial conditions w, 

F^[t<D<oo]< Ct-P/'^. 

Lemma 10. There is a di > such that, 

[V <oo]< l-5i. 

Proof. Without loss of generality, r = 0. Now, take the worst case scenario where w 
has M particles at each site x < — 1. By Lemma [HI 

Fu,[V = oo] > e-2A^/"'(i _ e"^'Y^U^^^{l - e^"^-)^'^ = 5i > 0. 

□ 

Lemma 11. Suppose that M > 4. There is a S2 > such that for all initial 
conditions w with at least a — 1 particles at the rightmost visited site r 

P„, [U <oo]<l- 62. 
Proof. We can also assume that r = 0. To estimate Fw[U = 00] below, note that 

F4U = 00] = p[n^^,{EU < ^/«'}]- (43) 

Let n G N and < e < l/(2a') and define G to be the event that each random 
walk Y^x,i) with < x < n, moves M + 1 steps to the right before time e. When G 
happens, Uk < l/(2a') and hence Yl'j=i^j < k/{2a') for all fc < n + [M/aJ := n' . 
Note that G n J^^'^^ Vj < k/a'} dGDH where 

H := nr=n'+i {E '^j < k/a' - n'e} (44) 

Furthermore, G and H are independent so Fw[U = 00] > P[G]P[H]. Now 

00 k 

pm < E ^[ E ^ ^/"' - ^'^] 

k=n'+l j=n'+l 

and letting = Vj - E[uj] = i/j - 1/a, for 2 < p < M/2, with /? := 1/a' - 1/a > 

-=^'+1 > k/a' - n'e] < i?[(E U'+i " n'e)-P 

<GkP/^{kp-n'e)-P, (45) 

where in the last inequality we used the same estimates explained between displays 
(ESI) and |(Sni). Taking e = /?/(2n') gives kP/'^{kf3 -n'e)-P < G'k'P'"^. Hence < 
CEfe^rt'+i SO as long as p > 2 (which is possible since M > 4) we obtain that 

P[H^] < 1 — (52 < 1 for sufficiently large n. Choose such an n < 00, and note that 
for this n, P[G] > > as well. □ 
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Lemma 12. Suppose that M > 4. There is a 5 > such that for all initial conditions 
w with at least a — 1 particles at the rightmost visited site r 

F^[D <oo]<l-6. (46) 

Proof. Since U and V are independent by part (iii) of Lemma 01 P«;[-D < oo] = 
1-F.^[U = oo]¥4V = oo]. 

□ 

Lemma 13. There is a C < oo such that for every initial condition w with r = 0, 
and any M' > M, 

[rt > M't] < Cexp{-Ct}. (47) 

Proof. Note that rj is a process on Z, increasing by one whenever a particle jumps 
to the right from r. Since there at most M particles there, the maximum jump rate 
is M. The Lemma then follows from standard estimates on Poisson processes. □ 

Lemma 14. For each p < M/2 there is a C < oo such that 
Proof. Let us first write, 

oo 

KSo [f^l > t| f/ = CX)] = ^Pa5o [Sk>t,K = k\U = CX^]. 

k=l 

Applying recursively the strong Markov property to the stopping times {Sj : j > 1} 
we see that for every > 1, 

Faso [Sk >t,K = k\U = (x,]<{l- 5f-\ 

where 5 > is given by Lemma IT^ For any ^ > we therefore have, 

FaSo [/^i > t| f/ = oo] < Ei=i^a5o [t<Sk<oo\U = oo] + 6~\l - 6)'. (49) 

Let 1 > 7 > and consider the event 

Ak = {roi - < t^, rD2 -rs^ <t^,..., rD^-^ - rs^_-, < t^} 

On Ak we have rs,^ < k{L + t^). Since < r^, if [/ = oo, then rt > [a't\ for all 
t > 0. Therefore, on Akn{U = oo}, 

[a'Sk] < k{L + t^). (50) 

Hence for t > (£(L + 1^) + l)/a' and k < £, 

FaSo [t< Sk<oo,Ak\U = oo]=0 

and therefore 

Kso [t < Sk < c^\U = oo]< Faso [Al, 5fc < oo| [/ = oo] . (51) 

By Lemma [TT| since FaSo [?7 = oo] > ^2 > 0. So for some C < oo, the right hand side 
of is bounded above by 

C Eti ^aSo [ro. -rs,>t"',Sk<oo]. (52) 
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Now let M' > M and 

Pa5o [rD. - rs, > t\ Sk < cx)] =¥aSo [ra. - rs, > 5fc < cx), A - S^ < f/M'] 
+Kso [rn. - rs. > t\ Sk < oo, A - > tyW] 
< Kso h.+tyM' - rs. > t^] + Fa5o [tVM' < A - Si < oo] . (53) 

Note that in the last equation we used the fact that Sk < oo implies that Di — St < oo 
for i < k. By the strong Markov property and Lemma [TSl 

Pa5o K+n/Af - rs, > t^] < Cexp-Ct^ 

By the strong Markov property and Corollary for each p < M/2, 

P,5o [tVM' < Di - Si <oo]< Ct"'<P/\ 

Choosing i = Clogt with C = p(21og(l - (5)"^)"^ from ^ and the previous 
estimates we obtain (EHl). □ 



Corollary 2. Let ki be the first regeneration time of the stochastic combustion pro- 
cess. 

a) For M > A, ¥.aSo [^i| U = oo] < oo, and ^aSo ['"kJ U = oo] < oo. 

b) For M > 8, EaS^ [i^i\U = oo] < oo, and EaS^ [rl^ \U = oo] < oo. 

Proof. The statements for ki follow from Lemma El and those for r^j from Lemmas 
Uniandlll □ 

6. Limit theorems 

In this section we use the renewal structure to prove Theorem the law of large 
numbers and the central limit theorem for r^. The argument for the law of large 
numbers follows that of Sznitman and Zerner in ^Hl; developed in the context of 
multi-dimensional transient random walks in random environments. The argument 
for the central limit theorem is from |17| . 

6.L Law of Large Numbers. We consider the stochastic combustion process 
started with an initial condition r = 0, ??(0,0) G {!,..., M} and r]{0,x) G 
{0, . . . , M}. We will prove that a.s. 

1- n Ea5o[r Kj\U = oo] 

lim — = V := ^ '\ ,' f. (54) 

t^^ t EaS,[i^i\U = oo] ^ ' 

Let us first note that by Proposition |S1 we have a.s., 

lim — = EaSo[Ki\U = oo], and lim — = EaSo[ri^i\U = oo]. (55) 

n— »oo n n—>oc n 

For t > 0, define nt := sup{n > : Kn < t}, with the convention that kq = 0. Note 
that (ESI ensures that nt < oo a.s. and by definition we also have, Km < t < k^^+i, 
and liuit^oo i^nt = c«. It follows from this inequality and the limit (|55)l that, a.s., 
limt^oo nt/t = l/EaSoinilU = oo] and hence almost surely, 

hm r^ /t = lim (r«,„ ) (k„, /t)=v, (56) 

I— >oo t— >oo ' 
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Now 

lim < lim ''^"*+^ ~ = 0. (57) 

t—^00 t t^OO t 

from This proves the law of large numbers. 

6.2. Central limit theorem. Starting with the same initial conditions we consider 

Bt := e'/' (r,-it - e-\t) , t > 0, 

Define 

-Rj := J^Kj+i - ''kj - {i^j+i - i^j)v, j > 0, 

and denote for m > the partial sums := J2]=i ^j- 
For any < t < T < oo, 

\Bl - e^/^Sn^^J < 2e^/^ sup (r-^^+i - r^J + 2ue^/^ sup (k„+i - k„). 

0<n<n^^„ij,j 0<n<n^^_ij,j 

For every u > we have, by Proposition [HI 

IP'a5o[suPo<n<n^^_lyj e^/^(K„+l - K„) > u] 

< P^^Jki > +n-2e(te-i + 1)E,5„ [Kfl(Ki > e-i/2^)|f/ = oo], 

which by part (b) of Corollary a.s. converges to as e — > 0. Hence, in probability 

sup e^/^(K„+i - Kn) 0. (58) 

0<n<?i|^^_l yj 



and similarly 



sup e^^^{rK„+, -r^J 



Hence, Bf — e^/^E„ converges to in probability, uniformly on compact sets of 
t. From Donsker's invariance principle, we know that y^S./^ converges in law to 
a Brownian motion with variance E5o[(rKj — kiv)'^\U = oo], where Ss,s > 0, now 
stands for the linear interpolation of Y!,m,m > 0. From the previous proof we have 
liuit^oo nt / 1 = l/KaSo[Ki\U = oo]. Since e~^kf^-i is increasing in t, the convergence 
is uniform on compact sets of t. This, together with the convergence in law of e^^^S./^ 
to a Brownian motion with variance Ea5g[(rKi — kiv)^], implies that e^/^S^ is 
tight in the Skorohod topology, and that its finite dimensional distributions converge 
to the finite dimensional distribution of a Brownian motion with variance. 



2 ._ ^aSoiirKi - KiVy\U = oo] 



E,5j/^i|t/ = oo] ■ ^^''^ 
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6.3. Non-degeneracy of the variance. It suffices to prove that, for some a! < 

13 <v, 

IPa5o Vk, = L, < I [/ = oo] > 0. 

Now, 

Pa5o = L, <Ki,U = oo]> FaSo [Lr^ <Si<U,Does,= oo] 
The right hand side we can write as 

Ea5o[l{Lp~' <Si< U)KaS,[l{V o 63, = oo)l([/ o 63, = oo) I 
Given ^5^, V o and U o 63-^ are independent, so 

EaSoiUV o 03, = 00)1([/ o 03, = 00) I .FsJ] 

= Pa5o [1^ o = 00 I J-sJlPa^o [C/ o = OO | ^3,]]- (60) 

But since by Lemmas [TI71 and [TTI we have Pa5o[^ ° ^S'l = 0° I -^^i]] = Pa^olt^ = 00] > 
62 > and Pa5o[y o 83-^ = 00 I .F^J] > 61 > we have, 

P„5o[L/3-i < 5i < o = 00] > <5i<52Pa5o[i/3~' <Si<U]. 

But it is easy to check that FaSo[L(3~^ < Si < U] > 0. In fact, it is enough to lower 
bound this probabihty by the probabihty that one of the random walks at site 
moves to site L in a time t such that L(3~^ < t < L{a')~^, and then stays at site L 
between time t and time L{a')~^, while all other random walks between sites and 
L do not move at all during the time interval [0, L{a')~^]. 
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